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We give a systematic method to calculate higher derivative corrections to low-energy 
effective theories of solitons, which are in general nonlinear sigma models on the moduli 
spaces of the solitons. By applying it to the effective theory of a single BPS non-Abelian 
vortex in U{N) gauge theory with N fundamental Higgs fields, we obtain four derivative 
corrections to the effective sigma model on the moduli space C x CP^~^. We compare them 
with the Nambu-Goto action and the Faddeev-Skyrme model. We also show that Yang- 
Mills instantons/monopoles trapped inside a non-Abelian vortex membrane/string are not 
modified in the presence of higher derivative terms. 



Email addresses: meto(at)riken . jp, toshiaki . f uj imori (at)pi . inf n. it , 
iiitta(at)phys-h.keio .ac.jp, ohashi (at) gauge . scphys . kyoto-u .ac.jp, sakai (at) lab . twcu .ac.jp 



Contents 

1 Introduction 

2 The non-Abelian vortices 

3 Preliminary: a particle in R" 

4 The effective theory 

4.1 Derivative expansion 

4.2 General formulas for the effective Lagrangian 

4.2.1 The zeroth order effective Lagrangian 

4.2.2 The second order effective Lagrangian 

4.2.3 The fourth order effective Lagrangian 

5 Correction to the single vortex effective Lagrangian 

5.1 Comparison with other models 

5.2 Instantons trapped inside a non-Abelian vortex 

5.3 Higher derivative terms and massive modes 

6 Mass deformations 

6.1 Higher order corrections from mass terms 

6.2 Kink monopoles 

6.3 Q-solitons 

7 Summary and Discussion 
A 0(3) model 

B Gauged linear sigma model and higher derivative corrections 

C Regularity of the second order solution 



1 Introduction 



Solitons are smooth localized solutions of nonlinear partial differential equations, and they are 
ubiquitous in diverse fields in physics and mathematics. Topological solitons are smooth lo- 
calized solutions of field equations in quantum field theories, and their stability is ensured by 
topological charges. A particularly interesting class of topological solitons is Bogomol'nyi-Prasad- 
Sommerfield (BPS) solitons [1], which saturate minimum energy bounds called the Bogomol'nyi 
bounds. They are the most stable solutions with given topological numbers and satisfy first order 
partial differential equations called BPS equations. Prominent examples of such BPS solitons 
are instantons in Euclidean Yang-Mills theory [2], 't Hooft-Polyakov monopoles [3] in the Bo- 
gomol'nyi limit [1], and Abrikosov-Nielsen-Olesen (ANO) vortices |1] at the critical coupling (a 
standard reference is [5]). Recent examples of BPS topological solitons in gauge theories are BPS 
kinks (domain walls) in Abelian [Hj and non-Abelian [7] gauge theories. 

While they preserve half of supersymmetry and are quantum mechanically stable in super- 
symmetric gauge theories [5] on one hand, one of the most important features of BPS solitons is 
that there exist no static forces among them and consequently a continuous family of configura- 
tions with degenerate energy is allowed. As a result, generic solutions of BPS solitons contain 
collective coordinates, called moduli parameters as integration constants. The moduli parame- 
ters of BPS solitons parametrize the space of whole solutions of BPS equations, called the moduli 
space. Although there is no static force among BPS solitons, they non-trivially scatter each other 
when they are moving. While such dynamics was studied by computer simulations, it is very 
difficult to study it by means of a fully analytic approach. The seminal idea of Manton [9] is that 
when solitons move slowly, their dynamics can be described as geodesies on their moduli space. 
This low-energy approximation is now called the moduli, geodesic, or Manton approximation. 
After his work, the moduli spaces are recognized as the most important objects associated with 
BPS solitons [5]. The Atiyah-Hitchin metric on two BPS monopoles is one of the most prominent 
examples [TU]. By examining geodesies on it, one can study the scattering of two BPS monopoles 
and find some interesting phenomena such as the right angle scattering in head on collisions. The 
moduli space of BPS monopoles was further studied in [TTl [12]. The moduli space and soliton 
scattering were studied for ANO vortices in the BPS hmit [HI [H [151 [161 HH [TU [19] and for BPS 
domain walls [20] • More recently, the moduli space dynamics has been successfully applied even 
to 1/4 BPS composite solitons such as domain wall networks [2T] and strings stretched between 
parallel domain walls [22] . Although there has been a lot of remarkable progress in study of the 
moduli space of BPS solitons, the analysis of their dynamics has been restricted to the leading 
order of the low-energy (small- velocity) limit, except for a few examples [231 121] • 

Non-Abelian vortices were found recently in A/" = 2 supersymmetric U{N) QCD [23 I2B] . 
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Although there has been a lot of progress in the study of non-Abelian vortices [271 EH] , here we 
concentrate on their moduli space. A single non-Abelian vortex has moduli space C x CP^^^ 
where the former denotes the position and the latter, called orientational moduli, comes from 
the SU (N) vacuum symmetry broken in the presence of the vortex [251 ES] • One of the most 
important facts about the orientational moduli is that confined monopoles and trapped Yang- 
Mills instantons in the Higgs phase are realized, respectively, as kinks [29] and lumps [30l [31] 
in the vortex world-sheet CP^~^ sigma model. This world-sheet effective theory provides a 
physical explanation for the relationship between BPS spectra in two-dimensional A/" = (2, 2) 
sigma models and four-dimensional M = 2 supersymmetric QCD [521 E3]- The moduli space 
of multiple vortices with full moduli parameters was completely determined (without metric) in 
[3ll[28l[35]; The moduh space for k separated vortices is a /c-symmetric product (C x CP^"^)^ /Sk 
of the single vortex moduli space [M] while the whole space is regular. General formula for the 
moduli space metric and its Kahler potential were given in [SB]- The metric of the moduli 
space on the coincident vortices was found [371 138] and was applied to low-energy dynamics of 
two non-Abelian vortices in head-on collision, reconnections of two non-Abelian cosmic strings 
j37j and flux matching of vortex- monopole composite [SH]. Recently the moduli space metric of 
multiple non-Abelian vortices has been finally obtained for well separated vortices [ID], and their 
low-energy dynamics has been examined [4T] . 

The purpose of this paper is to propose a systematic method to study dynamics of BPS 
solitons moving with higher velocities beyond the Manton approximation; we give a general 
scheme to calculate derivative corrections and apply it to obtain four- derivative corrections to 
the low-energy effective theory on a single non-Abelian vortex. In the spirit of the low-energy 
effective action, the leading order terms can be obtained by integrating out massive modes. The 
lowest order terms take the form of nonlinear sigma models with two derivative terms for light or 
massless fields, typically Nambu-Goldstone modes. A famous example is the chiral Lagrangian or 
chiral perturbation theory for pions which are the Nambu-Goldstone bosons of chiral symmetry 
breaking. The next-leading terms consist of four or higher derivative corrections, which are the 
focus of our study. In order to obtain four derivative terms, we solve equations for massive modes 
and eliminate them order by order. As a concrete example, we consider a single non-Abelian 
vortex. We obtain four-derivative corrections for the translational moduli C, the orientational 
moduli (the CP^~^ model) and their mixing terms. 

In the literature a different expansion is known: an expansion from the Nambu-Goto action 
[42j in powers of (the inverse of) the width of strings. The effective action of a string in the thin 
limit can be described by the Nambu-Goto action [13]. The finite- width correction in general 
takes the form of extrinsic curvature squared, which is called stiffness or rigidity of string [H]. 
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Dynamics of a string with the correction term was studied [15]. The exphcit calculation of the 
rigidity term was performed in the Abelian Higgs model by many authors in the context of cosmic 
strings and QCD strings [16]. This expansion is in powers of #(5) — where we denote 

the number of derivatives and fields as #(9) and i^{X), respectively. The leading term, the 
Nambu-Goto action, is the zeroth order term in the series, and thus contains the same number of 
derivatives and fields. The next-leading term, the stiffness term, contains two more derivatives 
than the number of fields starting from One advantage of this width expansion is 

that even the leading term, the Nambu-Goto action, contains infinite number of derivatives so 
that it can describe large fluctuations or bends of a string. However the expansion in powers of 
i^{d) — i^{X) is not applicable for dynamics of moduli in general; for instance multiple solitons 
or even a single non- Abelian soliton with internal degrees of freedom. — On the other hand, our 
expansion is purely a derivative expansion commonly used in the literature of field theories. It is 
valid at low energy and even the leading term, the nonlinear sigma model, can describe dynamics 
of multiple solitons. 

In the present paper, we will discuss the higher derivative corrections to the low-energy 
effective theory on a non-Abelian vortex world-volume. By a symmetry argument, one can 
easily write down the generic form of the higher derivative terms and find that the corrections 
cause instability for the lumps (sigma model instantons) which are identified with the Yang- 
Mills instantons trapped inside a vortex. According to the Derrick's scaling argument, the size 
of a lump expands in the presence of the generic higher derivative terms. Since the lumps are 
responsible for non-perturbative effects in the vortex world-sheet effective theory, one may think 
that the correspondence between the BPS spectra in 2d and 4d [321 [33] would be modified by 
the higher derivative corrections. However, as we will see in Sec. l5.2[ the lump solutions are still 
stable if the higher derivative terms have a specific form. In this paper, we will calculate the 
explicit form of the higher derivative terms and find that they do not modify the lump solutions. 

This paper is organized as follows. We first review non-Abelian vortices and their moduli 
space in Sec. [21 After illustrating our method for the derivative expansion in a simple example 
of classical mechanics in Sec. [31 we develop a systematic method to obtain derivative corrections 
to the effective action for non-Abelian vortices in Sec.[4l In Sec. 15. II we compare our result for 
the effective action of a single vortex with previously known models, the Nambu-Goto action 
and the Faddeev-Skyrme model [17]. In Sec. 15. 21 we discuss that the four derivative terms do not 
modify the Yang- Mills instanton solutions trapped inside a non-Abelian vortex [301 [31] • In Sec. [HI 
we discuss higher order corrections to the vortex effective action in a mass deformed model and 
show that the classical BPS spectrum is not modified by the higher order corrections. Sec. [7] is 
devoted to summary and discussion. 
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2 The non-Abelian vortices 



In this section, we briefly review the non-Abehan vortices and summarize the basic tools to 
describe the moduli space of BPS configurations. 

Let us consider the U{N) gauge theory in {d+ l)-dimensional spacetime with gauge field PF^ 
and A'f = Higgs fields H {N-hj-N matrix) in the fundamental representation. The Lagrangian 
of our model takes the form 

^2 



Tr 



(2.1) 



where g is the gauge coupling constant and is the Fayet-Iliopoulos parameter. We use the al- 
most minus metric t]^,, = (+1, —1, —1, ■ ■ ■ , —1). Although we can choose different gauge coupling 
constants for the U{1) and SU{N) parts of the gauge group, we set them equal for notational 
simplicity. Our notation for the covariant derivative and the field strength is 



(2.2) 



As is well known, the Lagrangian can be embedded into a supersymmetric theory with eight 
supercharges. The vacuum condition is solved by 



H = vl 



N- 



(2.3) 



This vacuum expectation value (VEV) of the Higgs field completely breaks the gauge symmetry, 
whereas the following SU{N)c^p diagonal symmetry remains unbroken 



H UcHUf, ul = UFe SU{N)c+F- 



(2.4) 



Let us consider the non-Abelian vortices in this color-flavor locked Higgs vacuum. We assume 
that the vortices are localized in the z = Xd-i + ixd plane. The topological sectors of the field 
configurations are classified by the vorticity, i.e. the magnetic flux of the overall U{1) gauge 
group 



k 



2tt 



dzAdz TtF, 



keZ. 



(2.5) 



For a given vorticity k, the tension (the energy per unit volume) of static configurations is 
bounded from below as 



> 2TTV^k. 



AVy,{V,^HH^) + 2vHF, 



(2.6) 
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where we have assumed that ViH — at the infinity \z\ — oo so that the third term does not 
contribute to the energy. This Bogomol'nyi bound is saturated if the following equations are 
satisfied 

V,H = Q, iF,, = ^iv'lN-HH^). (2.7) 

These equations are the BPS equations for non-Abelian vortices and describe the configurations 
of static BPS vortices which are minimal energy configurations in a given topological sector. Since 
the same equations can be derived by imposing the condition that one half of supersymmetry 
is preserved, Eq. (12. 7p is also called 1/2 BPS equations. We can show that any of the BPS 
configurations is a solution of the following full equations of motion of the system 

= V^V^H + ^{HH^ -vHn)H, (2.8) 
= ^V^'F^, + i[H{V,H)^ -{V,H)H^]. (2.9) 

In order to describe the moduli space of the BPS configurations, it is convenient to write the 
BPS solution as 

H = v S~^Ho, Wg = -zS^^dS, (2.10) 

where Hq and 5* are A^-by-A^ matrices. Then, the first BPS equation in Eq. (12. 7p becomes 

d,Ho = 0. (2.11) 

Therefore all the entries of the A^-by-A^ matrix Hq are arbitrary holomorphic polynomials. The 
matrix function S G GL{N,C) is determined from the second equation in Eq. (12. 7p . which can 
be rewritten into the following equation for Q = SS'^ : [HI [3i 

4 



g2y2 



dsino.n-^) = HoH^n-^ - 1^. (2.12) 



For a given Hq, this equation can be viewed as a non-linear differential equation for Q. The 
boundary condition for Q should be chosen so that the energy density vanishes at the spatial 
infinity \z\ oo. In general, this condition is solved by 

n — y HqHI (2.13) 

The matrix S can be determined from the solution Q uniquely up to the unphysical U{N) gauge 
transformation S — )■ Su{x) (with u{x) G U{N)). Solving Eq. (12.120 for an arbitrarily chosen 
Hq{z), we can obtain a BPS vortex solutions through Eq. (I2.10p . In other words, the matrix 
Ho[z) classifies all the BPS configurations. Hence the holomorphic matrix Hq[z) is called "the 
moduli matrix" and Eq. (I2.12p is called "the master equation" for vortices. 
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The matrices Ho{z) and Q are invariant under the original U{N) gauge transformation. On 
the other hand, the original fields H and are invariant under "the ^^-transformations" defined 
by 

Ho{z) ^V{z)Ho{z), n^V{z)nV{z)\ V{z) eGL{NX)- (2.14) 

This is the gauge symmetry of the master equation Eq. (12.121) which does not change the phys- 
ical quantities. Therefore, the solutions of the BPS equations are completely classified by the 
equivalence classes of the moduli matrix 

Ho{z) ~ V{z)Ho{z). (2.15) 

This implies that the parameters contained in an appropriately gauge fixed moduli matrix Hq 
can be interpreted as the moduli parameters of the BPS configurations. 

Let us next discuss how the moduli matrix Hq determines the topological sector of the cor- 
responding solution. From Eq. fl2.10p . we can see that the magnetic fiux of the overall U{1) is 
given in terms of Q by 

Tr F,-, = -idsd, log det fi. (2.16) 

By using the Stokes' theorem and the boundary condition Eq. fl2.13p . we find that the vorticity 
Eq. (12. 5 p is given by 

k = {dzd,- dzd-,) log det{HoHl), (2.17) 

where S^^ denotes the clock- wise circle at infinity \z\ — oo. This equation shows that the order 
of the polynomial det Ho{z) corresponds to the vorticity. 

As an example, let us consider the case of single vortex configurations. The simplest single 
vortex solution can be obtained by embedding the single Abelian vortex solution in the upper-left 
corner of the A^-by-A^ matrices H and 

H = [ , W, = [ ' . 2.18 

The other single vortex solutions can be obtained by acting on the embedded solution with the 
SU{N)c+F symmetry 

H U^HU, Ws U^W.U, U e SU{N)c+F- (2.19) 

Since the vortex embedded into the upper-left corner breaks SU{N)c+f SU{N — 1) x f/(l), 
the internal vortex moduli space is the complex projective space 

CP^-i ^ SU{N) 

~ su{N-i)xuiiy ^^-^^^ 
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By modeling out the unbroken SU{N — 1) x U{1), we can fix the unitary matrix U as 

\ l;v-i J \ Y-2 ) \ h ) 

where the parameters h = {b^, ■ ■ ■ ,&^~^) can be interpreted as the inhomogeneous coordinates 
of CP^~^ and X and Y are given b}|l] 

X = l + P-b, Y = lN^i + b®bl (2.22) 

The physical meaning of the orientational moduli can be seen from the magnetic flux 

F,, = X ( \ ] . (2.23) 

1 + V -b b®b^ J ^ ^ 

— * 

Thus the orientational moduli b determines the U{1) generator in which the vortex magnetic flux 
is embedded. 

The moduli matrix Hq and the solution of the master equation Q corresponding to the single 
vortex configurations are given by 

Ho{z) = Vi ]U, Q = Vi ]V\ (2.24) 

where ip is a real profile function satisfying 

^ dAi^ = l- \z- Zl'^e"^, (2.25) 



g2y2 



with the boundary condition — )■ log |z — Zp. The complex parameter Z in the moduli matrix 
Hq can be interpreted as the position moduli of the vortex. Although the holomorphic matrix 
V{z) can be an arbitrary element of GL{N,C), it will be convenient to use the matrix of the 
form 

Viz) = i 1 • 2.26 

V Y2 J \ J 

This ^^-transformation is chosen so that the moduli matrix Hq{z) takes the form 

Hoiz) = f ° 1. (2.27) 



b 1 



The important point is that the ^^-transformation is completely fixed so that the moduli matrix 
is holomorphic not only in z but also in the moduli parameters Z and b. In general, there exists 



^The square root of the matrix Y is given by 1^=*=^ = (1 - P) + (1 + \b\'^)^i P with P = |6|"26 (g) . 



such a fixed form of the moduh matrix in each coordinate patch of the moduh space. The V- 
transformation between a pair of the fixed moduh matrices induces a coordinate transformation 
from one patch to another. For instance, in the case of = 2, there are two fixed forms 





^0 = ^ . I , ^0=1. . I • (2.28) 
These two matrices are related by 

= VHo, V={ \ " 1 , 6' = f (2.29) 




Thus, the induced coordinate transformation is the standard transition map between the in- 
homogeneous coordinates of CP^. As this example shows, the coordinate fransformations are 
holomorphic, so that the moduh space is a complex manifold. 

For a general winding number k, we can fix generic moduli matrices to the following form 



Ho{z)=\ ^ ^ ^ , ei{z) = [[ y y (2.30) 

y bieiiz) + b2e2{z)--- + bkek{z) l^^i J j^i^ 

The parameters Zj and bj {I = 1, ■ ■ ■ , k) are position and orientational moduli of I-th vortex 
which cover a local coordinate patch of the moduli space of the vortices Mk- As in the case of the 
single vortex, the moduli matrix Ho{z) is holomorphic with respect to the moduli parameters. 
This fact will be important when we derive general formulas for the second and fourth order 
effective Lagrangian in section HI 



3 Preliminary: a particle in 

Before studying higher derivative corrections to the vortex effective action in field theories, we 
first explain our basic strategy in a simple system of a particle in M" with the following Lagrangian 

L = -i.i-V(x), x(^) = (x^x^...,a;"). (3.1) 

In a minimum energy configuration, the particle stays at the bottom of the potential V where 
the gradient of V vanishes 

grad\/ = 0. (3.2) 
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Now let us assume that the potential V has flat directions which are parameterized by 0*. Then 
we can define "the moduli space" of the minimum energy configurations by 



M = {x.^'^\(j)') \ giadV = } , (3.3) 

and 0* can be interpreted as "the moduli parameters". Since one can shift the particle to any 
points on the bottom of the potential without loss of energy, we can assume that the particle 
slowly moves along the moduli space Ai for sufficiently small excitation energy. This motion of 
the particle can be represented by the moduli parameters 0*(t) which weakly depend on the time 
t, that is 

x(t)=x(°)(0^(t)). (3.4) 

Then, the low energy effective Lagrangian of the particle can be obtained by substituting Eq. (13.41) 
into the original Lagrangian Eq. (13.11) 

Le. = + = -l-(x(°)) + ^x(°) . x(0), (3.5) 
where = —V{x.''^^) is the constant value of the potential at the bottom. The second order 

(2) 

effective Lagrangian L^g can be rewritten by using the metric gij on the moduli space Ai as 

4? = J9^,^'^■ (3.6) 
The metric is given by the inner products of the basis of "the zero modes" 

(7., = Ax(o). (3.7) 

The effective equations of motion for the moduli parameters take the form of the geodesic equation 

4>' + Ty'<P' = 0, (3.8) 

where the connection is given by 

1 d 

r}fc = 2^'^ i9ij,k + gikj - gjk,i) = /$;-^$a:. (3.9) 



This approximation is valid if the velocity of the particle is sufficiently small so that the time 
derivative is much less than the typical mass scale of the massive modes determined from the 
Hessian matrix H of the potential V defined by 



02 

rHfx^^hi = — V 

L V J jab dx^dx^ 



(3.10) 

AO) 
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Next, let us consider higher derivative corrections to the effective Lagrangian by talcing the 
massive modes into account. To this end, we first introduce a correction to (I3.4p by adding small 
fluctuations to "the slowly moving background" x*^°)(0*(t)) as 

x(t) = xW(</.Xt))+5x(^). (3.11) 

Since the motion of the particle along the flat direction is already represented by the moduli 
parameters 4>\t), we impose the following condition to avoid the double-counting of the degrees 
of freedom of the zero modes: 

5x-$i = 0. (3.12) 

This means that the fluctuations 5x is orthogonal to the tangent space of Ai, namely the fluctu- 
ation 5x contains only massive modes. Then, the original Lagrangian Eq. (13. ip can be rewritten 
as 

L = l(5+l2 + 5L + A^(5x-$,), (3.13) 
SL = mi(°)-5i+y5x-(5x-^(5xH(5x+--- , (3.14) 

where we have introduced the Lagrange multipliers A* to impose the constraint Eq. f l3.12p . Note 
that there is no linear term in the Taylor expansion of the potential since the background satisfles 
grad V^(x'^'')) = 0. Now let us consider the expansion of 5x with respect to the time derivative dt 

5x = x(2)+x(^) + --- , x(")~0(9r). (3.15) 

Note that 5x can have only terms with even numbers of the time derivatives due to the reflection 
symmetry t — )■ —t. There is no zeroth order term since the fluctuation 5x vanishes for the static 
conflgurations. Correspondingly, the Lagrangian 6L can also be expanded as 

(5L = L(^)+L(6) + ... . (3.16) 

The fourth order Lagrangian contains the background x*^"-' and the second order fluctuation x*^^) 

L(^)=mi(°)-i(2)_lx(2)Hx(2). (3.17) 

Note that x*^^^ does not contribute to the fourth order Lagrangian since gradl^(x(o)) = 0. The 
massive modes x^^^ can be eliminated by solving their equation of motion 

Hx(2) +mi(°) = A*$i. (3.18) 

To determine the Lagrange multiplier A*, let us take the inner products of the both hand sides 
and the zero modes $j 

m$j-x(°) = gijX\ (3.19) 
10 



where we have used the fact that the zero modes satisfy 

= A[grady(x(°))] = 0. (3.20) 

We can show that the Lagrange multipher is proportional to the second order equation of motion 
Eq.dMD 



X = m g'^ $j ■ ±^^> = m (^0* + r^;(/.>' j . (3.21) 

Substituting back A* into Eq. fl3.18p . we obtain the following equations of motion for the massive 
modes 

Hx(2) = -mPi(°\ (3.22) 
where P is the projection matrix which project out the zero modes 

Px(o) ^ ^(0) -^^i(i{o) . $.)$^, (3.23) 

Since the right hand side of Eq. f l3.22p does not contain the zero mode directions, the matrix H 
can be "inverted" as 

x(2) = _mGi(°), (3.24) 

where G is the matrix satisfying 

HG = P, GP = PG = G, G$, = 0. (3.25) 

Substituting the solution for the massive modes Eq. f l3.24p into Eq. (I3.17p . we obtain the following 
fourth order effective Lagrangian 



2 2 ^ ^ ^ ^ \ d(j)' dcj)'' J 



(3.26) 



Let us see a simple example of a particle in M^. We assume that the potential is rotationally 
symmetric V = V^(|x|) and has a minimum at |x| = tq, namely V'{ro) = and V"{ro) > 0. The 
static configurations are parameterized by the moduli parameter 6 as 

,n^ / rn COS 9 \ , , 

xW = V ) . (3.27) 
Y ro sm u J 

The second order effective Lagrangian is given by 

l2 = = -rle'. (3.28) 
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The corresponding equation of motion for the moduh parameter is ^ = and describes the 
particle rotating around the circle at |x| = tq. The matrices H and G are respectively given by 



H 



V'\r,) P, G = y"(ro)-ip, 



(3.29) 



where P is the projection operator 

cos^ 9 sin 9 cos 9 
sin 9 cos 9 sin^ 9 

From Eq. (13.261) . we obtain the following fourth order Lagrangian 



(3.30) 



-(4) 
^eff 



(3.31) 



Even if we take into account this fourth order Lagrangian, the equation of motion is not modified 
^ = and solved hy 9 = ut + 9q. On the other hand, the relation between the angular velocity 
a; and the angular momentum / is modified as 



/ = mrl [l + 2ma;V"(ro)~^] uj. 
This is just because the rotation radius is increased by the centrifugal force as 



To 



ro + mro 9"^ ^"('^o)"^- 



(3.32) 



(3.33) 



The shift of the rotation radius can be also seen in the solution for the massive mode 



,{o) 



m 



cos 9 
sin 9 



(3.34) 



Therefore, the higher derivative term gives the correction from the massive mode which is slightly 
shifted by the motion of the zero mode (see Fig.[T]). 




Fig. 1: The shift of the rotation radius. 
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4 The effective theory 



4.1 Derivative expansion 

In this section, we discuss the effective Lagrangian for non-Abelian vortices by generahzing the 
method of the derivative expansion discussed in the previous section. To deal with the vortex 
positions and orientations Z/, 6/ (/ = 1, ■ ■ ■ , A;) on an equal footing, we combine them into a 
set of complex moduli parameters (z = 1, ■ ■ ■ , kN = dime A^fc), and assume that the moduli 
matrix Hq{z) is always holomorphic in cff like Eq. f l2.30p . 

The zero mode fluctuations along the vortex world- volume are described by the vortex moduli 
(ff promoted to fields which weakly depend on the world-volume coordinates 

(a = 0,1,--- ,rf-2). (4.1) 

These moduli fields induce fluctuations of massive modes around the weakly fluctuating vortex 
background 

= + (4.2) 

W,{x^) = lyf (z,z,0(x")) + 5W,{x^), (4.3) 
Wa^x^") = + SW^i.x^'), (4.4) 

where H^^'^ and H^i"^ are the BPS vortex background Eq. fl2.10p depending on the world-volume 
coordinates through the moduli fields 0(x°). We assume that the excitation energy of the 
fluctuations are much less than the typical mass scale of the massive modes. Then we can expand 
the induced fields with respect to the derivative by assuming that 

9„ < gv. (4.5) 

Note that gv is the unique mass scale controlling the mass of the bulk fields, inverse width of 
the vortex and the mass scale of the massive modes localized on the vortex world-volume. The 
induced fluctuations are expanded with respect to the derivative da as 

5H = + ^(4) ^ . . . ^ ^4_g) 

= Wf^ + Wf + ■ ■ ■ , (4.7) 
= Wi^^ + Wi^^ + -- - . (4.8) 

These fluctuations can be determined by solving the equations of motion (12. 8p and (12.91) order- 
by-order. Note that odd (even) order equations of motion for 6H and 6W2 {^Wa) are trivial 
due to the reflection symmetry x°' — )• —x°'. Then, the effective Lagrangian can be obtained by 
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eliminating the induced fluctuations from original Lagrangian (12. ip expanded with respect to the 



derivative dn 



C = £(0) + £(2) + £(4) + . . 



(4.9) 



As we will see, the Lagrangian is quadratic in the induced fluctuations up to the fourth order in 
the derivative da-, so that we can use the linearized equations of motion to eliminate the induced 
fluctuations. In order to determine the fluctuations, we have to specify the boundary conditions 
for them. In the singular gauge, the vortex background takes the following form at the infinity 

\z\ — )■ oo 



■gv\z\'' 



(0) 



0{z 



(4.10) 



Since this vortex background is independent of the zero mode fluctuations 0(x°) at infinity, the 
minimal excitations induced by the zero modes should vanish at infinity. Therefore, we impose 
the following boundary conditions for the induced fluctuations 



5H — > 0, 



5W„, 



0, 



(in the singular gauge). 



(4.11) 



4.2 General formulas for the effective Lagrangian 



4.2.1 The zeroth order effective Lagrangian 



The zeroth order Lagrangian contains only the background fields 



£(0) ^ 



9^ 



N 



(4.12) 



Substituting the solution fl2.10p and using the boundary condition Eq. f l2.13p . we find that the 
zeroth order term of the low-energy effective Lagrangian is given by the sum of the tension 
of the individual vortices 



(0) 
eff 



T 



2'KV' 



(4.13) 



4.2.2 The second order effective Lagrangian 



The second order Lagrangian takes the form [251 ISSj (see also [331 SHI EOl [3l]). 



£(2) = Tr 
where F^^^ cind VaH'^^^ are given by 



(4.14) 



daH^^^ + tWi^^H^^\ (4.15) 
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Note that the terms containing the second order fields if*-^^ and W^i^^ do not contribute to the 
second order Lagrangian since these terms are proportional to the background equations of motion 







Tr 



6 

6H 



d 



'd+l 



X 



£(0) 



BPS background 



(4.16) 



The dynamical degrees of freedom in the second order Lagrangian are not only the zero modes 
0*(x") contained in the background fields {H^^\ W^i^"*) but also the first order fiuctuations Wa^\ 
We can eliminate Wa'^ by using the equations of motioij^ 



n \ ^ zn '^z^ za 



(4.17) 



By using the fact that the moduli matrix Hq{z) is holomorphic in the moduli parameters 0* and 
the matrix Q is the solution of the master equation Eq. f l2.12p . we can check that the solution is 
given by 



where we have defined the differential operators Sa and 5^ by 

Sa = da 



d<p^ 



(4.18) 



(4.19) 



The solution Eq. (14.18^ satisfies the boundary condition Eq. ( 14. lip since the asymptotic form 
of the matrix 5* in the singular gauge is S* — )■ ifo(^;0*) and the moduli matrix Hq{z, </)'■) is 
holomorphic with respect to 0*. Substituting the solution ( I4.18P into the second order Lagrangian 
(I4.14p . we obtain the following formal expression of the second order effective Lagrangian 



(2) 
eff 



c/^x 5lTr 



This effective Lagrangian gives a natural Kahler metric on the moduli space of vortices 



(fx --^Tr 



(4.20) 



(4.21) 



Note that the above formulas for the first order solution fl4.18p and the moduli space metric are 
not invariant under the generic V^-transformation Eq. (12.141) . However they are invariant under 
the V^-transformation respecting the holomorphy of Hq{z) as we pointed out at the end of Sec.[21 

In order to obtain the explicit form of the moduli space metric, one needs to know the 
solution of the master equation Q. However, no analytic solution of ^2 has been known even 
for the minimal winding vortex, and hence it is quite difficult to obtain the explicit form of 



^ Note that the orthogonahty condition between zero modes and fluctuations is trivial for wi^^ since Wa = 
in the vortex background. 
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the effective Lagrangian in general. Nevertfieless, we can obtain an exact form of the second 
order effective Lagrangian for the single vortex since the degrees of freedom Z, b in this case are 
nothing but Nambu-Goldstone zero modes [30]. By substituting the solution Eq. fl2.24p and using 
the boundary condition for the profile function ip — )■ log \z — Zl"^, we can show that 



C 



(2) 
cfT 



9 



(4.22) 



where q^-^ is the Fubini- Study metric on CP^ ^ 



9^ 



FS 



(4.23) 



Another example of analytic moduli space metric for higher winding vortices has been obtained 
for well-separated vortices [iO] . 

Note the moduli space metric can be rewritten as inner products of zero modes as in the case 
of the example discussed in the previous section (see Eq. (13. 7p ). Before closing this subsection, 
let us review the zero modes and the derivative operator defining them [25]. The physical zero 
modes in a BPS background are defined as the solutions of the linearized BPS equations for the 
fluctuations around the background. With a suitable gauge fixing condition, the linearized BPS 
equations can be written as 



A<l> 



3-H 



(0)t 




(4.24) 



where the subscript r denotes the fact that H^^^ acts as right multiplication and P| and "D" are 
covariant derivatives with the zeroth order gauge field which act on the fundamental and adjoint 
fields, respectively. There exists one zero mode $j for each moduli parameter 0* and the set 
{$j} (z = 1, ■ ■ ■ , dime A^) forms a basis of the zero modejf|. In terms of 1), S and Hq, the basis 
can be written as 



\ 9 



s 



(4.25) 



^ The linearized equations of motion in a static BPS background can be summarized into the form of A^^A<I> = 0. 
Here, the Hcrmitian conjugate of the operator A 



turns out to have no zero mode 



A.. I 

At$ = 



a fjm 



afjio) 



<f> = 0, so that gives a full set of the zero modes. 
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By using the solution for wji'^ given in Eq. (14.18^ . these zero modes can be summarized as 



2f(1) 



(4.26) 



Then, the moduh space metric Eq. f l4.2ip can be rewritten as inner products of the physical zero 
modes 



(4.27) 



where we have defined the hermitian inner product for pairs of fields in the fundamental and 
adjoint representations 



(4.28) 



4.2.3 The fourth order effective Lagrangian 



Let us now calculate the fourth order effective Lagrangian by solving the linearized equations of 
motion and eliminating the massive modes. The fourth order Lagrangian takes the form 



£(4) ^ 



9^ 



4 



I),-i/(2)+,H/(2)i7(o) 



(4.29) 



Here, terms proportional to Wa'^\H^^^ and , wi^^ automatically vanish due to the same mecha- 
nism as Eq.f l4.16]) . It will be convenient to combine H'^^^ and W^'^^ into a column vector 



$(2) = 



ij(2) 



Note that the fluctuations transform under the gauge transformations as 

^(2) _^ ^(2) ^ ^^^(0) 



(4.30) 



(4.31) 



where A is an arbitrary hermitian matrix of order d^. In order to eliminate this unphysical 
degrees of freedom, let us impose the following gauge fixing conditions for the the fluctuations 

.2 

(4.32) 



^p^l^f + ^i7(2)^(o)t = 



This constraint is equivalent to the condition that the second order fluctuation is orthogonal 
to the unphysical gauge zero modes $a 







($A,$(2)) + (/l.C.), 



a ^ 



(4.33) 
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Then, the terms containing {H^'^\ W^i^"*) in the fourth order effective Lagrangian are sum- 
marized as 



4( A<l>(2) A<l>(2) 



(4.34) 



where we have introduced the Lagrange multipher A* to impose the condition that the second 
order fluctuation is orthogonal to the physical zero modes $j, as we have done in Eq. (13.131) 
in order to separate the massive modes from the zero modes. The linearized equations of motion 
for and l^i^^ can be written as (cf. H ^ A^A in Eq. fIXTSD ) 



4A^A$(2) +r'„(9"0^$i) = X'^i. 



(4.35) 



Let us first determine the Lagrange multipliers A*. By taking the inner products of the zero 
modes $J and the both hand sides of the linearized equation Eq. (14.351) . we find that 



(4.36) 



where we have used 



k 



, AtA$(2)^ = ^A<l>j- , A$(2)^ = 0. (4.37) 



We find that A* = is nothing but the equation of motion for 0* with the second order Lagrangian 
Eq. (l4.20p . Then, the linearized equation Eq. (I4.35P becomes 



AtA<l>(=^) + -P 
4 



0, 



where P the following projection operator which projects out the zero modes 



(4.38) 



(4.39) 



It seems that in order to solve Eq. (I4.38P in terms of ^^^^ , we need to know the explicit form of 
the background BPS solution . However, we can find the following formal expression for A$(^) 
without solving the BPS equations 



A$(2) 



2 

2^ 



gv 



d 



V 



7Tt-l 



\ 



S 



(4.40) 



where Vj is the covariant derivative on the moduli space, which acts on (9n~^f2) as 



d d 
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(4.41) 



We can confirm that Eq. fl4.40p satisfies the hnearized equation Eq. f l4.38p by checking that the 
following equivalent equation is satisfied 

1 



4 



A 



(AAt)A$(2) 

where the operator AA"!" takes the form 

-I)/p/ + ^(i7(o)tij(o)) 



AA^ 



2! 
4 









(4.42) 



(4.43) 



The solution A^^^^^ is unique since the operator AA''" obviously has no zero mode and hence it is 
invertible. Although it has apparent singularities due to the factor Hq~^, we can show that the 
solution Eq. f l4.40p is smooth everywhere (see Appendix [C] for the proof). 

It is in general not easy to solve Eq. fl4.40p in terms of However, the explicit form of 

A^*^^-* is sufficient for the purpose of determining the fourth order Lagrangian since it can be 
rewritten as 

1 



(4) 
efF 



d xTi 



29 



2 ^^al3 ) 



(4.44) 



By using the first order solution fl4.18p . the first term in the fourth order Lagrangian fl4.44p can 
be calculated as 



Tr 



1 



2^7 



2 



9^ 



5f^(l]5^]fi-i)5t["(f]5/31f]-i) 



(4.45) 



Substituting f l4.45p and fl4.40p into (I4.44p . we obtain the following form of the fourth order 
effective Lagrangian for the moduli fields 



,{4) 

-eff 



where the tensors A^^u and B^.u are given by 



^ijkl 



— I cPxTi 
9^ 



— I d xTi 

9^ 



9 9 

gZyZ 



(4.46) 

(4.47) 
(4.48) 



5 Correction to the single vortex effective Lagrangian 



Applying the formula Eq. fl4.44p to the single vortex solution Eq. (12.240 . we obtain the following 
fourth order effective Lagrangian for a single non-Abelian vortex 

T_ „_„,o 47r 



gi^Zd^^Z - Ki''^{d^Zd^Z) + c{9lfd''h^d^h^ 



■ (5.1) 
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where T = 2tiv'^ is the vortex tension and c is a constant given by 



— / d^x il-\z- Z^e-'^f ~ 0.830707 x 

ZTT J ' g^v^ 



(5.2) 



In summary, we have obtained the following derivative expansion of the effective Lagrangian for 
a single non-Abelian vortex 



eff 



+ - 



+ 



T 
Air 

7 

Atic 



i-ldc^zd'^z-hdaZd'^z 

z o 



{gfd^b'diV) 



(5.3) 



5.1 Comparison with other models 

Let us compare our result, Eq. 05.31) . of the four derivative terms in the effective action for a 
single non-Abelian vortex with those in different context. First we compare the translation 
modes Z with the Nambu-Goto action, and second we compare the orientational modes 6* with 
the Faddeev-Skyrme model. 



Translational zero modes and mixed terms 

First, by setting the orientational moduli 6* to zero in Eq. ( 15. 3p . we have the effective action 
for the translational modes Z as 



d^-^x 1 - -d^Zd'^Z 



1 



(5.4) 



which is precisely the one of an Abelian (ANO) vortex. At this order, this effective Lagrangian 
coincides with the Nambu-Goto action [12] given by 



^1 



NG 



-T j d^-^x y/-det(-7„/3), 



where 7q,/3 is the induced metric on the world-volume given by 



Vap - - [daZdpZ + daZdiiZ) . 



(5.5) 



(5.6) 



In Eq. (15.31) . there are the four derivative terms containing both the translational zero modes Z 
and the orientational zero modes U. These are precisely the terms appearing in the derivative 
expansion of the following action describing the <CP^~^ sigma model on the vortex world- volume 

An 



S, 



eff 



d'^~'xJ-det{--fap) 



-T + -rP{g^fd^UdpV) + 0{dl] 



(5.7) 
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Orientational zero modes 

Next, we consider the orientational zero modes in the internal space. For simplicity, let us 
restrict ourselves to the case of = 2 in which the internal moduli space is CP^. By setting 
Z = in Eq. fl5.3p . we obtain the following effective Lagrangian of the orientational zero modes 
with the higher derivative correction 



An 

7 



(5.8) 



Now let us compare this with the Skyrme-Faddeev model which is a CP^ model with a four- 
derivative term [17]. To this end, we formulate the CP^ model (at the leading order) by two 
complex fields h = (/i^, h"^) charged under U{1) gauge symmetry. By introducing auxiliary gauge 
field tta and scalar field A as Lagrange multipliers, it can be written as 

C = (da + iaa)h(da — iaa)h^ — \ ( hh^ — ^] (5.9) 

V r / 



The variation of A gives the constraint hh^ = which can be solved by h = \/ ^'^J\t\i (1? ^)- 



Then the variation of gives 

i bdab — bdjj 



\b\' 



(5.10) 



Substituting these back into the original Lagrangian (15. 9p . the CP^ model (at the leading order) 
is recovered. 

The field strength of the gauge field (I5.10p is 

tap = Oaap-dpaa = -I ^ i^p^^ . (5.11) 

Then the Faddeev-Skyrme term [17] can be written as the field strength squared aq^ 

. .aP , {dab(rbf-\dabcFb\'' 

Uf - 2 (i + |5|2)4 • (5-12) 

This term is quadratic in the time derivative and does not coincide with Eq. (15. 8p . The other 
term containing four time-derivatives appears in an A/" = 1 supersymmetric extension of the 
Skyrme-Faddeev term [51|: @ 

(dabd^'b)'^ 



\b\ 



2\4: ■ 



(5.13) 



This term is also called the baby Skyrme term in d = 2 + 1 [SO] • 
^ This term also arises when one constructs the Faddcev-Skyrme-like model as the low-energy effective theory 

of pure SU{2) Yang-Mills theory [52], and hence is called the Gies term in that context. 
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As shown in Appendix[B], this higher derivative term can be obtained by adding a higher derivative 
term to the original Lagrangian f l5.9p . In general the fourth order terms are summarized as 

.(4) _ , {d^bd'^br _ {2c, + C2){dJd-br-2c,\d^bd'^b\' 

^general " CiJq,/?/ + '^2 ^-j^ |^|2^4 (^1 _|_ |5|2)4 " K^-'-^) 

This reduces to the four derivative term in Eq. ( 15. Sp when 

C2 = -2ci = c. (5.15) 

The condition C2 = — 2ci is precisely the condition for A/" = 1 supersymmetry [SI]. This must be 
the case because there remains A/" = 1 supersymmetry (four supercharges) in the vortex effective 
theory, because vortices are 1/2 BPS states in supersymmetric theories with eight supercharges. 



5.2 Instantons trapped inside a non-Abelian vortex 

In this section, we consider d = 4 + 1 dimensions where vortices are membranes with 2+1 
dimensional world-volume. Besides the vortices, there also exist Yang-Mills instantons which are 
particle-like solitons in d = 4 + 1 dimensions. But they cannot exist stably in the Higgs phase 
which we are considering. Instead, they can stably exist inside the world- volume of a non-Abelian 
vortex [5DII5T] . Instanton- vortex composite configurations are 1/4 BPS states in supersymmetric 
gauge theories with 8 supercharges. In the supersymmetric CP^~^ model with 4 supercharges, 
at the leading order of the 1/2 BPS vortex effective theory, instantons can be regarded as 1/2 
BPS sigma model lumps [30] which are point-like solitons in (i = 2 + 1 dimensions. Here we show 
that instantons (= lumps) are not modified even if we include derivative corrections found in the 
previous sections. 

For simplicity, we set Z = and consider the case of = 2 (the CP^ sigma model) with the 
Lagrangian (15. 8 p for the orientational modes. In c? = 4 + 1 dimensions, the vortex is a membrane 
which has coordinates (x° = t, a;^,x^). We parametrize the two spatial coordinates by complex 
variables w = + ix"^ and w = x^ — ix"^. We discuss only static configurations in this model. 

First we recall the lump solution in the model without the four derivative correction. The 
energy density can be written as 

p.(2) _ 4vr f ^ d^bd'^b _ 8n f 2\d^b\' + {\d^b\' - \d^,b\') 



The first term is positive semi-definite > and the second term can be rewritten as 

Ani f d,„b d,T,b — d,nbd,„b , , Am f db A db 87r^ , 

dw Adw = — / -—— = -^k, (5.17) 



9' J (l + l&IT 9' J {l + m' 9' 
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where the integer k is the degree of the map from the vortex world-volume to the target space 
CP^. This is the topological charge for lumps: 7r2(CP^) = Z. The energy (15.161) is bound from 
below by this topological charge I, and the Bogomol'nyi bound is saturated by the BPS equation 
for lumps, given by 

= 0. (5.18) 

The BPS solutions b{w) are holomorphic in w and satisfy the static equation of motion. If we 
fix the boundary condition as 6 — i- oo (w — )• oo), the solutions with I = are given by k-th. 
order rational maps 

. aow'' + aiw''-^ ^ hafc 

o[w) = - — ^- r — : • (5.19) 

aotf + aiw'^ ^ + ■ ■ ■ + flfc-i 

This is the general solutions for the 1/2 BPS sigma model lumps. 

Next let us consider the effect of the four derivative term 

(4) _ _647rc f , \d^bd^b\' 

The total energy can be rewritten as 



1 -4c- 



■k. (5.21) 



Again, the first term is positive semi-definite as long as the correction term is sufficiently small and 
vanishes if the BPS equation (15.181) is satisfied. We thus have found that there is no contribution 
from the four derivative term to the BPS configurations and consequently the lump solutions are 
not modified. The energy of lumps ^k agrees with that of instantons in the Higgs phase, which 
can be found by rewriting the energy of the original bulk theory to the Bogomol'nyi form for 1/4 
BPS configurations. Therefore, it is natural to conjecture that all higher order terms vanish for 
1/2 BPS lump configurations. 

By considering 1/2 BPS lumps on the CP^~^ model, we have obtained 1/4 BPS composite 
states of instantons inside a vortex. On the other hand, there exist another 1/4 BPS composite 
states of intersecting vortex- membranes [301 EI]- Vortices of one kind have codimensions in z- 
plane and extend to w-plane, and those of another kind have codimensions in w-plane and extend 
to 2;-plane. This intersecting vortices can be constructed by considering holomorphic maps of the 
translational moduli Z{w), instead of the orientational moduli b^{w) considered in this section. 
Now let us show that, in general, the fourth derivative terms do not modify the 1/2 BPS states. 
The BPS equation can be found by rewriting the second order energy as 

^iff = J d''xAgfjd^(P'd^^ + J igf^dcj)' A d^K (5.22) 
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The second term gives the area of a two cycle of the moduh space on which the map 0* wraps. 
The first term is positive semi-definite and vanishes for any holomorphic maps satisfying 



<9,7 



0. 



(5.23) 



On the other hand, it follows from the general form of the higher derivative corrections Eq. fl4.46p 
that the fourth order energy vanishes for the holomorphic maps 



0. 



(5.24) 



Therefore, the higher derivative corrections do not modify the 1/2 BPS configurations. 

There is a possibility that the holomorphic maps are unstable if the fourth order energy is 
negative. However, this would be an artifact of the truncation at the fourth order. For example, 
the translational part of the energy can be rewritten as 



= T I d^x2{l-\d^Z\'^)\dij,Z\^ + T j ^ {dw A dw + dZ A dZ) 



(5.25) 



The second term is a topological term and the first term vanishes for a holomorphic map Z(w) 
but it is not positive semi-definite. Although a holomorphic map Z{w) appears to be an unstable 
solution if I^^Zp > 1, such a case is beyond the validity of our approximation. We can show 
that all the holomorphic maps are stable BPS solutions if we consider the full order action, i.e. 
the Nambu-Goto action, whose energy can be rewritten as 



NG 



T j d^x [V(l + |9^Z|2-|a^Z|2)2 + 4|9^Z|2 - (1 + \d^Z\^ - |9^Zn 
+ T j ^{dw Adw + dZ Adz) . 
The first term is positive semi-definite and can be expanded as 



T d^x2 



1 - \d^Z\^ + \d^Zf{\d^Zf + \d^Zf) + ■ ■ ■ \d^Z\^. 



(5.26) 



(5.27) 



As this example shows, we should take into account the full order corrections to prove the stability 
of the holomorphic maps. 



5.3 Higher derivative terms and massive modes 

To see the physical meaning of the higher derivative corrections for the internal orientation, let 
us consider the equation of motion with higher order corrections in the N = 2 case. Assuming 
that the orientational zero mode h is independent of the spatial world-volume coordinates, we 
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can show that the following configuration satisfies the equation of motion even in the presence 
of the higher derivative terms : 

b = exp{iut). (5.28) 

This solution corresponds to an excited state of the non-Abelian vortex whose orientation is 
rotating along the equator of CP^ (Fig.|2]). The rotation of the orientation induces one of three 



b = co 




b = 



Fig. 2: The rotation of the orientation 
components of the SU (2) c+f conserved charges, which is given by 

The relation between the angular velocity u and the conserved charge ps (shown in Fig.|3]-(a)) 
is modified by the higher derivative term, as in the case of the angular velocity and the angular 
momentum of the particle (see Eq. f l3.32p ). 

By numerically solving the full equations of motion Eqs. f l2.8p and f l2.9p . we can see that 
the size of the excited vortex is slightly larger than that of the static configuration (see Fig.|3]- 
(b)). This is analogous to the case of the particle discussed in Sec.|3l the higher derivative term 
corresponds to the correction from the massive mode (vortex size) which is slightly shifted by 
the conserved charge. 

6 Mass deformations 

6.1 Higher order corrections from mass terms 

In this section, we consider higher order corrections from supersymmetric mass deformations of 
the original theory. Without breaking the supersymmetry, we can deform the original model by 
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P3 



8 r 




(a) The relation between and u (b) The shift of the vortex size 



Fig. 3: (a) The relation between the angular velocity oj and the conserved charge pa for g = v = 1. Compared to 
the second order relation pa = p-w (dashed line), the fourth order relation p3 = ^-(1 + (solid line) shows 

the better agreement with the numerical result (dots) obtained by solving the full equations of motion Eqs. (|2.8p 
and (|2.9p . (b) The energy density distributions in the transverse plane (z-plane) for g = v = 1. The size of the 
excited vortex with ui = 0.5 (solid line) is slightly larger than that of the static vortex with w = (dashed line). 



adding the following terms to the Lagrangian Eq. (12.1 

1 „ . „ 1 



•^adjoint+mass Tr 



9 ^9 



(6.1) 



where S/ (/ = 1, ■ ■ ■ , ra) are the real adjoint scalar fields in the vector multiplets and M/ are 
mass matrices in the Cartan subalgebra of SU{N)p 

Mi = mj a = diag{mi^o, mj^i, ■ ■ ■ , mj^N-i)- (6.2) 

If all the masses are non-degenerate, the SU{N)p flavor symmetry is maximally broken to the 
Cartan subalgebra of SU{N)p. These mass deformations can be obtained by dimensional re- 
ductions from the six dimensional model on x (S*^)" with the following twisted boundary 
conditions around periodic dimensions 

Hix^", Oi + 2txRi) = H{x^, ej)e^^'^'^^'. (6.3) 

Ignoring the inflnite tower of the Kaluza-Klein modes, we obtain the mass deformation Eq. (16. ip 
with the following identiflcation 

Eiixf") = -We.ix^"), H{x^')e''^'^^' = Hixf^^Oj). (6.4) 

The mass terms do not change the VEV of H (given in Eq. ( 12. 3 p ) , while they induce those of the 
adjoint scalars 

(S/) = Mj. (6.5) 
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In this mass deformed model, the color-flavor global symmetry SU{N)c+f is explicitly broken to 
the Cartan subgroup f/(l)^~^. Hence the orientational moduli, which were the Nambu-Goldstone 
zero modes of SU{N)c^p, are lifted by a potential induced by the mass terms. The potential on 
the moduli space can be calculated by finding the minimum energy configuration and evaluating 
the energy for each values of the moduli parameters. This can be done perturbatively with 
respect to the masses mj^A. As a zeroth order configuration, we consider the vortex solution 
Eq. fl2.10p satisfying the BPS equations without the mass deformations. Then we can determine 
the corrections to S/ by solving their equations of motion 



2 

~i2 



(v^V^Jlj + [Sj[Sj, S,]]) = (SjiJ - HMi)H^ + H{H^T.i - MjH^). (6.6) 

Due to the fact that S/ originate from the higher dimensional gauge fields We^- , we can solve the 
equations of motion in the similar way to the gauge fields Wa 

= M + i{6iS^S^^^ - S^^d\S) + 0{m\ mdl), (6.7) 

where we have assumed that the derivatives da and the masses m are of the same order 

m ~ da- (6.8) 

The differential operators Sj and 6\ are defined by 

where k] are holomorphic Killing vectors on the moduli space which are the following linear 
combinations of the Killing vectors of the unbroken symmetry U{1)^~^ 

k\ = uiiC- (6.10) 

For example, the holomorphic Killing vectors on the moduli space of the single vortex solutions 
are given by 

k\ = i{mi i — mj o)b^ (no sum over i). (6-11) 

Inserting the solution T,^P into the deformation terms Eq. (16. ip . we obtain the effective po- 
tential of the form 

V}S^=g,jkW,. (6.12) 

This is the sum of the squared norm of the Killing vectors k] evaluated with respect to the moduli 
space metric g^j given in Eq. (14.211) . 
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The potential of this form can also be obtained by dimensional reductions from the four 
dimensional vortex world-volume action (the effective theory of the vortex in six spacetime di- 



mensions 



(fx gfj 



(6.13) 



As in the case of the bulk theory, we impose the twisted boundary condition for the moduli fields 

0^(x", 01 + 2tiRi) = 0^(x", ei)e'^^^'''l (6.14) 

Keeping only the lowest modes and evaluating the kinetic terms in the effective action, we obtain 
the effective action with the potential Eq. fl6.12p 



J, 



{4-n)d 



eff 



(6.15) 



Therefore, the procedure of the dimensional reductions and the calculation of the effective action 
are commutative. 

The higher order corrections to H and Wz can also be determined from the equations of 
motion with the mass deformations. They can also be solved in a similar way to the case without 
mass deformation as 



2^? 



gv 



\ 



H, 



S 



(6.16) 



Then the higher order corrections are obtained by substituting the solution into 



(4) 
eff 



d xTt 



(1) vWj2 



4{^A<I>(2) , A<l>(2)y (6.17) 



In the case of the effective Lagrangian of a single vortex, the induced terms in take the form 



(2)+(4) 
off mass 



Air 



+ 



9 

4ttc 
Attc 



9 



{gfl^jdJ)=){glfk^jd^V) + c.c. 
, {9fmMfkWj). 



(6.18) 



These corrections terms can also be obtained by using the dimensional reductions from the four- 
dimensional world-volume with the higher derivative corrections. 
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6.2 Kink monopoles 



In this section we consider a single 1/2 BPS non-Abelian vortex-string in four spacetime dimen- 
sions and its low-energy effective theory in two spacetime dimensions. Besides the vortex, an 
another important topological soliton, 1/2 BPS monopole, arises in four dimensions. However, 
since the gauge symmetry is completely broken in the Higgs phase, the monopole cannot stay 
alone and must be accompanied with the vortex-string. If the vortex-strings are attached to 
the monopoles from the both left and right hand sides, the composite vortex-monopole state 
is indeed a stable BPS state preserving 1/4 supersymmetry (two supercharges) [29]. Such con- 
fined monopoles can also be viewed as kinks in the two-dimensional vortex effective theory. For 
example, let us consider N = 2 case with only one mass term 



M 



- m, 
2^ ' 



-m 



(6.19) 



The orientational part of the second order effective Lagrangian of a single vortex-string in the 
mass deformed U (2) theory takes the form 

47r9«69"6-m2|6|2 



[a 



0,1 . 



(6.20) 



9' {i + m' ' 

In this sigma model, the orientational moduli are lifted by the potential and there are discrete 
vacua at 6 = 0, oo. In the (1 + 1) -dimensional vortex world-sheet, we can consider a BPS kink 
interpolating between the discrete vacua 

'Ait \dxb — mh\^ 



E. 



(2) 
eff 



dx 



a, 



9' {i + WY 

where a is the moment map of the unbroken U{1) symmetry defined by 

27rm 1 - 



9' l + W 



(6.21) 



(6.22) 



The BPS equation can be easily solved as 

dxh = mb, — )■ b{x) = exp [m(x — Xq) 



i9] 



(6.23) 



where Xq and 6 are kink position and phase moduli, respectively. The mass of the kink is 



Ml 



BPS 



a{x — oo) — a{x — )■ — oo) 



(6.24) 



As expected, this BPS mass of the kink coincides with that of the monopole in (3-1-1) dimensions. 

To see the effect of the higher order corrections to the BPS kinks, let us take into account 
the forth order corrections in the effective theory 



(2+4) _ ^ 

9 



2 



dabd^b-m'^\b\'^ \dj)d''b + m'^b 



2^212 



1^1 



2^2 



+ c- 



+ 



2\4: 



(6.25) 
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The energy of static configurations can be rewritten as 

Therefore, the BPS kink solution and its mass Eqs. (16.231) and (I6.24p are not modified by the 
higher order corrections. Thus, we again encountered the case where the topological soliton in 
the vortex effective theory is exact at the second order. This is consistent with the observation 
that the kink on the vortex corresponds to the monopoles attached by the vortices in the four 
dimensional full theory. 



'47r \dxb - mhy 



2^2 



\dj) + mby 



2\2 



+ 



(6.26) 



6.3 Q-solitons 

Let us again consider the effective theory of the single non-Abelian vortex in 4 + 1 dimensions. 
In the case of = 2, the effective theory is (2 + l)-dimensional massive CP^ sigma model. The 
effective Lagrangian at the second order is the same as Eq. (I6.20p with a = 0, 1,2. As is well 
known, from the Derrick's theorem, there are no stable lump solitons: the potential makes them 
collapse. However, one can still construct stable solitons by adding Noether charges to the lumps. 
It is the so-called Q-lumps [511 ESj . The Q-lump solutions have the same form as the CP^ but 
their phases are time-dependent. To see this, let us see the energy 

^(2) ^ ^ [d\ \^? + m'-m^h\' 



2 / 1 + / + (6.27) 



9' J {l + W? 

_ 47r /■ 2 \bTvrn^l±\du,b\''/2 

~ 9' J (1 + H 

with / and Q being the topological charge and the Noether charge associated with the U{1) 

global symmetry 

8ni f dbAdb Sn'^k 

(6.28) 



J (1 + |6|T ' 
The BPS equations are of the form 

b = ±imb, disb = 0. (6.30) 

Thus the solution is given by 

b{t,x\x^) = e^'"''b{w), (6.31) 

where b{w) can be any holomorphic function of w as in the case of the static configuration 
Eq. fl5.19p . Note that the Q-lump solution with the minimal winding number (6 = — ) has infinite 



30 



energy because its Noether charge density has an asymptotic tail ~ l/|A^w| [55]. The minimal 
configuration with finite energy, for example, is given hy b = — — . 

Let us next see the effects from the higher derivative corrections: Do the Q-lumps receive 
corrections from the higher order terms? Indeed, as one can easily check, the solution given 
in Eq. (16.301) solves the full equations of motion including the fourth order corrections. In the 
presence of the higher order corrections, the energy density can be rewritten as follows 



^^"""^ / <fx ^ . 1- -. (6.32) 



9^ J + 

Since is a holomorphic function in w, we can immediately conclude that E^^'^ = for the Q-lump 
solutions b = e*"^*6(w). Therefore, the BPS Q-lump solutions and their masses are not modified 
by the higher order corrections. 



7 Summary and Discussion 

We have proposed a systematic method to obtain higher derivative terms in the low-energy 
effective theories on solitons. We have applied our method to a single non-Abelian vortex and 
have obtained four- derivative terms in the C x CP^~^ model on the vortex world-volume. We 
have compared our four-derivative terms with the Nambu-Goto action and the Faddeev-Skyrme 
model. The action for the translational moduli Z coincides with the Nambu-Goto action and the 
terms for the orientational moduli 6* coincide with those for the supersymmetric extension of the 
Faddeev-Skyrme term. We have also shown that the contribution from the four-derivative terms 
disappears for 1/4 BPS states of instantons trapped inside a non-Abelian vortex and consequently 
the solutions are not modified in the presence of the four derivative terms. 

In this paper, we have derived the four derivative terms. In principle we can go on to any order 
in our formalism. The sixth order is considered to be prominently important. We have confirmed 
that the effective action of the translational zero modes Z at the fourth order is consistent with 
the Nambu-Goto action. On the other hand, in the width expansion from the Nambu-Goto 
action, the first correction term written as the extrinsic curvature squared starts from the sixth 
derivatives in the derivative expansion [H]. In field theory calculation, it seems that there is no 
agreement on the signature of that term [IB] • Therefore, we can in principle determine that term 
for the BPS case. 

In this paper we have studied local non-Abelian vortices which exist in the theory with the 
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number Nf of flavors equals to the number A^c of color. When the theory has more flavors, 
A'^p > Nq, vortices are called semi-local [56]. It is known in this case that the orientational zero 
modes of a single vortex is non-normalizable, i.e., the integration over the codimensions diverges 
j57] . unless the size modulus is zero and the vortex shrinks to a local vortex [SH]. For two vortices, 
the relative orientational zero modes are normalizable even with a non-zero size moduli, while 
the overall orientational zero modes are non-normalizable [58]. Although we can formally extend 
our method to semi-local vortices, we should check if there exists a divergence in four derivative 
terms even for normalizable moduli. Remember that in the derivative expansion, derivatives are 
assumed to be less than the lowest mass m of the mass spectrum in the vacuum. The existence 
of the vacuum moduli for A'^p > Nq implies that m = and hence the convergence radius of 
the derivative expansion seems to be zero. Therefore, it is interesting to see if there exists a 
special mechanism which justifies the derivative expansion for semi-local vortices. Non-Abelian 
vortices were extended to arbitrary gauge groups G in the form of {U{1) x G)/C{G) with the 
center C{G) of G [59]. Especially the cases of G = SO{N), USp{2N) have been studied in detail 
[60} [6T|. We can straightforwardly extend our analysis to the cases of arbitrary gauge groups but 
we should be careful to the normalizability since they are semi-local vortices in general. 

Although we have studied BPS vortices in supersymmetric gauge theories, our method to 
obtain higher derivative corrections by solving equations of motion for massive fields is robust 
and can be extended to non-supersymmetric theories. For instance, we can apply it to non- 
Abelian vortices in non-supersymmetric theory [62]. In reality, non-Abelian vortices exist in 
high-density QCD which may be realized in the core of neutron stars [63]. In this case, the low- 
energy world-sheet theory on the vortex is described by the bosonic CP^ model at the leading 
order [6l]. The four- derivative correction to it should be important especially for the fate of 
confined monopoles which have been recently shown to exist as kinks on the vortex [65]. This is 
because they are non-BPS and higher derivative terms do not vanish automatically, unlike BPS 
instantons discussed in Sec. 15.21 

Finally our method is general so that we can apply it to other BPS solitons such as domain 
walls, monopoles and instantons, or non-BPS solitons such as Skyrmions. In the same spirit, 
four derivative term in the form of the Skyrme term was obtained in the effective theory of 
non-Abelian domain walls [23], and four derivative terms for collective coordinates of a rotating 
Skyrmion were calculated [23]. 
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A 0(3) model 



From the isomorphism CP^ ^ ^ 0(3)/0(2), the CP^ model is equivalent to the 0(3) model. 
In order to see this equivalence, let us introduce a three vector n = (ni, ?t.2, ris) by 

n - I , ^ , ^ \a ( 



which satisfies the constraint 

= 1. (A.2) 

Conversely, /3 is the stereographic coordinate, given by 

/3 = !^1±^ = ^-"^ . (A.3) 
l + ns ni - in2 



The kinetic term becomes 



-d^n ■ d^n, (A.4) 



(1 + |/3|2)2 2 ^ 
and the field strength can be rewritten as 

ff,^ = n- {df,nx d^n). (A.5) 

Therefore the Skyrme-Faddeev term and the other four derivative term become 

f,ur = id.nxd^n)^ (A.6) 
-{d^n-d^nf, (A.7) 



(1 + I/3IT 4^^''* 



respectively. The total four derivative terms in Eq. f l5.14p can be rewritten as 

£4 = ci(<9^n X d^nY + ^{d,n ■ d^'nf. (A.8) 
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B Gauged linear sigma model and higher derivative cor- 
rections 



In this section, we consider a gauged linear sigma model with a higher derivative term which 
reproduces the orientational part of the higher derivative terms in Eq. (15. 3p . The model is de- 
scribed by charged scalar fields h = {hi, ■ ■ ■ , hiy) coupled to a U{1) auxiliary gauge field Oq. 
The Lagrangian is given by 

£ = V^hiV'^h)^ + dV^h{Vph)^V"h{V^h)^ + X{hh^ -v^) , (B.l) 

where d is a constant and A is a Lagrange multiplier for the constraint hh'^ = v"^, which is solved 
by 

h= JL^ ih,--- ,b^_„l). (B.2) 

The covariant derivative is defined by V^h = [da + iaa)h. The equation of motion for the 
auxiliary gauge field be solved as 

i hidj)i - dakbi 

This solution is independent of the parameter d. In other words, the higher derivative term in 
Eq. ( IB. II) does not change the solution for the auxiliary gauge field. We can easily show that 
Eq. fIB.Sp is the solution by using the following relations 

Vahh^ = h{V^hy = 0. (B.4) 

Substituting Eq. ( IB. 31) back into the Lagrangian Eq. ( IB. II) . we obtain 

£ = v^gfdJfd'^V + v^d{gfdj)'dpb>){glfd''h^d%^). (B.5) 

Therefore, ii v"^ = ^ and d = |^c, this Lagrangian coincides with the orientational part of the 
vortex effective action Eq. (15. 3p . 



C Regularity of the second order solution 

In this section, we show that there is no singularity in the first component of the solution 
given in Eq. (14. 40 p . Here we assume that Zj Zj (I J), so that the constant matrix 
= [{z - Zi)Hq^]^^^ is well-defined. Then the singular behavior around the J-th vortex 
position z = Zj can be written as 

+ 0(1). (C.l) 

z=Zi 
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Z — Zj 



d 



We can show that 



d 



J z=Zr 



dZi d 



dZ 



(C.2) 



- z=Zj 



where we have used the following identity for any function of the form f{z, z, (f)\ (jf) 
Then, the singular part can be rewritten into the following form 



4 z-Zi 



with Bi 



9 9 



7.1 



(C.3) 



(C.4) 



where we used the following relations which can be derived from the master equation Eq. ( I2.12p : 



2„,2 



TV 



vi/,, vp,a|[na,n-i] . =0. 



(C.5) 



\z = Zl 4 - \z = Zl 

Now let us use the explicit form of the generic moduli matrix Eq. fl2.30p . Then, is given by 



1 
-b 



(C.6) 



By assuming that the matrix Bi takes the form 

Pi {QiV 
^ ri sj 

we can rewrite the singular part as 



Bi 



(C.7) 



2,,2 



9 V 



4 z-Zj 



1 
-b 



d Pi- Zi {qi) 



n t 



ri 



Si - Zi 



We can show that Vi^{Zi — pi) = Vi-^qi = as follows. It has been shown that the matrix 
Bi is related to the moduli space metric as 

k 
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TTV 



nv 



7=1 

k 



dZi dZi 



Tr 



dH^^ dB^ 



(C.9) 



- z=Zj 



1=1 



dbi _ ^^T-, ir 



dZ, 



j=i 



where Dj^ = [dz^j{z)]^^^^. Since the vortex moduli space is a Kahler manifold, there exist a 
Kahler potential K, such that 

die 
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(C.IO) 



By comparing this equation with Eq.( ]C.9l) . we can read pi and g} as 

k 



d 



(Zi-Pi) 
dqi 



j=i 



db, 



—,\{{Zi-Zj)- 



die 



Since the moduh space metric is covariantly constant, it follows that 



0. 



(c.ii) 

(C.12) 



(C.13) 



This shows that there is no singularity in the solution Eq. fl4.40p . 
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